Exact convergence times for generation of random bipartite entanglement 



OO 
O 

o 

(N 
Oh' 

in 



Oh. 
a ■ 

^3 ; 



> 

p 

O 
OO 
O 



X 



Marko Znidaric 
Department of Physics, Faculty of Mathematics and Physics, 
University of Ljubljana, SI-1000 Ljubljana, Slovenia 

We calculate exact convergence times to reach random bipartite entanglement for various random 
protocols. The eigenproblem of a Markovian chain governing the process is mapped to a spin chain, 
thereby obtaining exact expression for the gap of the Markov chain for any number of qubits. For 
protocols coupling nearest neighbor qubits and CNOT gate the mapping goes to XYZ model while 
for U(4) gate it goes to an integrable XY model. For coupling between a random pair of qubits the 
mapping is to an integrable Lipkin-Meshkov-Glick model. In all cases the gap scales inversely with 
the number of qubits, thereby improving on a recent bound in [Phys. Rev. Lett. 98, 130502 (2007)]. 
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Introduction- Entanglement is a resource that can 
make quantum processes more powerful than classical 
ones. While a complete characterization of entangle- 
ment is a complicated task one might learn a great deal 
from properties of generic states, where by generic we 
mean states drawn randomly according to unitarily in- 
variant Haar measure, shortly random states. Random 
states have almost maximal bipartite entanglement [l[ 
and are a needed resource in some quantum information 
protocols like quantum dense coding 0] or remote state 
preparation In addition, they occur naturally dur- 
ing a time evolution of a sufficiently complex (quantum 
chaotic) system. Note though that detecting entangle- 
ment in random states could be difficult [1]. Central is- 
sue in quantum information processes is their efficiency. 
We want protocols that are faster than the best classical. 
Natural question therefore is, how difficult is it to gener- 
ate random quantum states? Knowledge about chaotic 
quantum systems might suggest that it should be possible 
to produce random states, as far as their bipartite entan- 
glement is concerned, in polynomial time, i.e., number of 
two qubit gates. On the other hand, one should be aware 
that to produce an arbitrary unitary transformation, and 
therefore a truly random state, an exponential number of 
two qubit gates is needed in general However, if our 
criterion is just to reproduce bipartite entanglement of 
typical random states, which is the case frequently stud- 
ied, less gates might be needed. 

Random protocol for generating random states is as 
follows: generate some pseudo-random sequence of two 
qubit gates applying it to an arbitrary separable initial 
state. After sufficient number of two qubit gates we will 
end up in a random state. Such random protocols have 
been numerically studied in [1, d, 0] ■ Recently [1| , ran- 
dom protocol for certain gates has been mapped to a 
Markov chain whose gap determines the convergence rate 
of a bipartite entanglement. It has been proved that the 
gap is lower bounded by ~ 1/n^ for CNOT gate and 
coupling between random qubits, rigorously establish- 
ing that only polynomial number of two qubit gates is 
needed to generate random bipartite entanglement. Nu- 



merical investigation [3] indicated that the bound is actu- 
ally ~ 1 /n, which was furthermore confirmed for random 
C/(4) gate and coupling between random qubits by an an- 
alytical study [§]. 

In the present paper we improve and extend on the 
bounds in ^ Q by analytically calculating the exact gap 
for C/(4), CNOT and XY gates and for various coupfing 
topologies like nearest neighbor with open and periodic 
boundary conditions and for coupling between all qubits. 
Note that analytical bounds for CNOT and XY gates 
are not known and that experimentally important case 
of nearest neighbor coupling has not been treated before. 
We find that in all cases the gap scales as '--^ 1/n, meaning 
that purity reaches level I ^ e after time scaling as ~ 
nlog (1/e). Method is easily generalizable to other gates. 

As a measure of bipartite entanglement we are going 
to use purity given by the trace of the square of the 
reduced density matrix pA[t) — trB[|V'(O)('0(*)|] ■ Ex- 
panding density matrix p — |-0(i))(?/'(t)| over products 
of local Pauli matrices, p — Cq cr" 
purity /(t) = trA[pi(i)] = 1^ Ea={a.o.} 4 
ct"' denotes Pauli matrix Ui S {0,x,y,z} acting on «th 
qubit, with the convention a'^ = 1. We want to calcu- 
late time dependence of purity for a protocol consisting 
of application of a random two-qubit gate Uij (t) acting 
on ith and jth qubits at step t, 1)) — Uij{t)\tlj{t)) . 

We shall consider two kinds of protocols: (i) Uij{t) is 
going to be a random U{A) gate, independent for each 
step, and (ii) Uij (t) is going to be a product of random 
single qubit gates V{t) and V'{t), independent for each 
qubit and at each step, and a fixed two qubit gate W, 
Uij{t) = Vi{t)V' {t)jWij. Averaging over random single 
qubit unitaries V and V from C/(2) one can arrive at 
the transformation law of coefficients Ca after one step 
of the protocol. If a two qubit gate W preserves prod- 
ucts of Pauli matrices, i.e., if W transforms a product 
of two Pauli matrices into a product of some other two 
Pauli matrices (apart from a sign), then the transforma- 
tion can actually be written for squares of Cq (for details 



, we get 
2 1+^^ Here 
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see the original derivation in [^), namely, 

c'{t + l)^McHt), M^i^Mf, (1) 

where by we denote a vector with components c^. As 
one can see, the protocol is described by a Markovian 
matrix M which is a sum of two-site Markovian matrices 
M^^-' between all coupled qubits. Two qubit matrix M'-^-' 
depends only on the two qubit gate W used in the proto- 
col. Its precise form will be given later. Markov matrix 
M has two eigenvalues equal to 1, corresponding to in- 
variant states being an identity operator and a uniform 
mixture of all basis states. We want to calculate the value 
of the 3rd eigenvalue, 1 — A, where A is the gap of Marko- 
vian matrix M. If this eigenvalue is nondegenerate purity 
will asymptotically decay as I(t) « exp (— iA) and will 
reach level I{t) — e after convergence time r = £• 
In the following we are going to calculate the gap of 
Markov chain M ([1]) for various 2-qubit gates, thereby 
obtaining decay time of purity. We will be in particu- 
lar interested in one-dimensional chains of n qubits with 
open (L — n — 1) or periodic (L = n) boundary con- 
ditions or with coupling between all pairs of qubits. For 
pedagogical reasons we shall give a detailed derivation for 
random U(4) gate with other gates being very similar. 

U(4) gate- Let us first have a look at the most sym- 
metric case when two-qubit transformations Uij{t)aTe 
independent random U(4) matrices. In this case 0, 0| 
two-qubit matrix M^"^^ is equal to Pigxie defined as. 
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F is here (m — 1) x (m— 1) matrix. We defined P to have 
an arbitrary dimension m because we are actually going 
to solve a more general situation with a two-site Markov 
matrix M*-^'' given by to x m dimensional matrix Pmxm- 
While dimension m can be arbitrary, we are in particular 
interested in cases when Pmxm acts on a tensor product 
of two local bases of dimension fc, that is to = fc^. If 
we are dealing with qudits with local operator basis of 
dimension (P we have m = d^, for instance, Piexie cor- 
responds to qubits and is the one we are most interested 
in. We are first going to study Piex 16 f^nd at the end just 
state results for some other dimensions. 

Any matrix can be decomposed into a sum of tensor 
products with the minimal number of terms, in quan- 
tum information sometimes referred to as the operator 
Schmidt decomposition [loj . 
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Such decomposition can be obtained by singular value de- 
composition where Kj are positive singular values while 



Aj , Bj are "columns" of unitary transformations occur- 
ring in singular value decomposition and act only on a 
single site. For Piexie such decomposition has only 4 
terms and due to the symmetry Bj can be chosen to 
be equal to Aj. Crucial observation is that all 4 x 4 
matrices Aj have the same two dimensional kernel (null- 
space) spanned by vectors (0, — 2/\/6, l/-\/6, l/\/6) and 
(0,0,-l/\/2,l/y2). Therefore, using unitary transfor- 
mation, Aj = UAjW, with 



U 




(4) 



A'j have a block form with nonzero elements only in the 
upper left 2x2 corner. Nontrivial local space is therefore 
only of dimension 2! Resumming these 2x2 blocks in Aj 
and Bj ([3]) one gets 
is exactly equal to XY hamiltonian 



a reduced P[q'^ ig of size 4x4 which 
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i6xl6 = ^(l®l + /lXY), 
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Spectrum of Markov chain M ^ will therefore be a 
union of eigenenergies of the full n site XY chain ([5]) and 
eigenenergies of all sub-chains obtained by dropping some 
sites and the corresponding couplings M^j connecting 
these sites with the rest. These sub-chain eigenenergies 
come from cases when an operator on the corresponding 
site comes from the kernel of Aj . 

For nearest-neighbor coupling we can use exact results 
for XY chain in magnetic field ll| to get the largest three 
eigenvalues of M for Pi6xi6- Although solutons are well 
known, for the sake of completeness and to properly treat 
boundary conditions we give details in the Appendix. For 
periodic boundary conditions the resulting gap is 



A = 



2(1 - ft. cos {■K/n)) 



(7) 



with h = 4/5. For open boundary conditions eigenvalue 
1 — A is doubly degenerate with the gap 



A = 



1 — /icos (j^/n) 



n- 1 



(8) 



In both cases the largest three eigenvalues come from the 
full chain of n sites as the eigenvalues corresponding to 
sub-chains are strictly smaller. Note that parameters of 
XY model satisfy + h? = \ for which the ground state 
becomes doubly degenerate [12|, corresponding to two 



invariant states of our Markov chain. Equations [7] and [8] 
provide exact expressions for the gap of Markov chain for 
periodic or open boundary conditions for any n. 
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Next, let us proceed with the case when the two-qubit 
gate is allowed between an arbitrary pair of qubits, that 
is when the sum in M ([1]) runs over all L — n{n — l)/2 
distinct pairs. Calculation of the gap, at least in the 
limit n oo, should be easier for such infinite range 
coupling as one could, for instance, use mean field ap- 
proximation applicable for infinite dimensional systems. 
We can actually do better though. First, one notes 
that for a general XYZ coupling, ft-xYZ = Jx'^f'^j + 
Jycrf crj + Jzcri<7j + + <^j), between all pairs of spins, 
^rcd ^ Y.i<ji(^'^ O 1 + hxYz)/L, we can rewrite this 
Markov matrix in terms of operators of a total spin, 

= UTn=l<^i)^ " =X,y,Z, 



prcd 
81x81 ~ 2 



2h 



— ( Jx'5'x + JySy + JzS^^ + 

d- ' : 1 1, (9) 



n{n — 1) 
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which is the Lipkin-Meshkov-Glick model (LMG) [l3|, 
solvable by a rather complicated algebraic Bethe 
ansatz ij]. Note that the dimensionality of the eigen- 
value problem has been reduced from ~ 4" for the orig- 
inal M to 2S + 1 ^ n for each total spin sector S. As 
we are mainly interested in the scaling of the gap we 
are going to calculate A in the limit of large n. Be- 
cause largest eigenvalues of A/''"* ([9]) come from the sec- 
tor of maximal spin S = n/2 we only have to calcu- 
late the difference between the two largest eigenvalues of 
LMG in this sector. There are various possibilities, per- 
haps the simplest one being by replacing quantum spin 
with a classical one, parameterized by a canonical pair 
ip and fi = COS0 as Sz = Sfj,, = S^/l -— cos tp and 
Sy = S'yl — /?sin iy9, expanding the resulting classical 
hamiltonian around its maximum to the lowest order in 
^ and ip, and subsequently quantizing the resulting har- 
monic oscillator hamiltonian, see e.g. [3]. We will just 
state the final result for the gap, which is for our param- 
eters of C/(4) gate (0) equal to 



on 



(10) 



Eigenvalue 1 — A is nondegenerate, while the next one 
(whose distance to the 2nd one is ~ 1/n?), coming from 
a sector with S = n/2 — 1, is n — 1 times degenerate. 
This degeneracy is simply due to ri — 1 time multiplicity 
of sector with total spin 5 = n/2 — 1 [3]. 

In a similar way one can treat also other dimensions 
of Pmxm- For instance, if local operator basis is of 
dimension k — 2, that is for matrix P4X4, one can 
again show equivalence with XY model hxy ®j this 
time with parameters 7 = 1/3 and h = 2\/2/3. For 
Pgx9 we have the equivalence Pg'^g = ^{1 1 + hxv) 
with 7=1/2 and h = ■\/3/2. Physically interesting 
case is also that of qutrits, for which the equivalence is 



+ hxy), with parameters 7 = 4/5 and 
h = 3/5. In all these cases the gap is given by eqs. (I7I8I) . 

CNOT gate- If W is CNOT gate two-site transforma- 
tion matrix Af(2) is equal to 0,1" 



DPa 



4x4, 



(11) 



where D is a, permutation matrix giving 
transformations of products of Pauli matri- 
ces by CNOT gate. Its matrix elements are 

Dt'+ic\t+4c = litr[WcNOTCr*cr'=W^(l^oyf7*'a'='] I = 

St'+4c'.f{t+4c) with Kronecker 6, /(0, 1, . . . , 15) = 
(0,1,14,15,5,4,11,10,9,8,7,6,12,13,2,3) and c, t 
denoting target and control qubits, respectively. We 
proceed along the same lines as for U(4) case. For 
CNOT we now have only three terms in the sum ([3]) for 
M'-^\ Kernel of matrices Aj and Bj is two dimensional 
and spanned by the same two vectors as for Pigxie- 
Therefore, rotation by unitary U, Eq. (jl]), will bring 
matrices Aj and Bj to a block form. Because CNOT 
lacks the full symmetry of U(4) there are now nonzero 
elements also in the lower left 2x2 block, in addition to 
an upper left 2x2 block. However, as we are interested 
only in right eigenvectors and eigenvalues we can again 
take only a sum of products of upper left 2x2 blocks, 
which is this time equal to the XYZ model. 



M. 



(2) 



1 



1 



7=l,/i=-. (12) 



XY gate- XY gate is given by Wxy|01) = 
-i|10),iyxY|10) = -i|01),WxY|00) = |00),iyxY|ll) = 
1 11), and was shown in [7] to be faster than CNOT 
or U(4) for nearest-neighbor couplings. The form 
of the two-site matrix M^^^^ is the same as for 
CNOT ()lip with the permutation matrix now given 
by A'+4c'.t+4c = St'+4c',g(t+4c) with .g(0, 1, . . . , 15) = 
(0, 11, 7, 12, 14, 5, 9, 2, 13, 6, 10, 1, 3, 8, 4, 15). Writing 
M*^^) as a tensor sum ^ we have 6 terms with the ker- 
nel of each matrix being again at least 2 dimensional and 
containing the same two vectors as for Piexw- Rotation 
by U (eq. SI results in matrices of the same structure as 
for CNOT gate, with the sum of products of nontrivial 
2x2 blocks being equal to XYZ model. 



18" 



1 , 2 

1, 7 = -, ft, = -. 
' ' 2' 3 
(13) 

For both CNOT and XY gate we obtain equivalence 
with the XYZ model in magnetic field. For nearest- 
neighbor coupling the model is not exactly solvable. Nev- 
ertheless, due to the existence of two invariant states of 
M we immediately know that the ground state of a fer- 
omagnetic chain is doubly degenerate for these param- 
eters. Chains also have a nonzero energy gap, meaning 
that the gap A will scale as ~ 1/n. For the case of 
coupling between all pairs of qubits the model is again 
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equivalent to LMG model ([9]). Using the same procedure 
as for U (4) gate we arrive at the gap 



A = 



4 

3n 



(14) 



which is the same for both CNOT and XY gate 7] . Inter- 
estingly, the gap is larger for CNOT and XY gates than 
for random ?7(4) gate (fTO|) . reflecting the fact that many 
random U{A) gates are only weakly entangling. Regard- 
ing degeneracies, for both gates the 2nd and 3rd largest 
eigenvalues (coming from S — n/2) are nondegenerate, 
while the 4th one (coming from 5* = ?i/2 — 1) is n — 1 
times degenerate for the same reason as for U (4) . This 
explains a seeming cutoff- like behavior observed in 8]. 

Conclusion- We have analytically calculated the gap 
of Markovian chain governing convergence of bipartite 
entanglement to that of random states. All calculations 
proceed by mapping transition matrix to various spin 
models. For C/(4) gate and nearest-neighbor coupling we 
obtained exact expressions for arbitrary n. For coupling 
between all pairs of qubits asymptotically exact expres- 
sion is obtained for U{4:), CNOT and XY gate. In all 
cases gap scales as ~ 1/n. The used method could be 
also employed for other gates or for qudits. 

Appendix- We want to find eigenenergies for XY 
model, H = T^, ®i with periodic or open bound- 

ary conditions [ll|. Note that for open boundary con- 
ditions magnetic field has strength h/2 on 1st and last 
spins. This is important as otherwise a doubly de- 
generate ground state for /i^ -I- 7^ — 1 splits into an 
exponentially close doublet. Using standard Wigner- 
Jordan transformation into fermionic operators Cj — 

■ ■ ■ a^_-f^{a^ — iaf)/2, one gets for periodic bound- 



ary conditions in spin variables, CTn+i ~ en, 

- h.c.) - 

Last term in 



En— 1 / 



H = 
hn + 



(_l)E'fe=i i+c,c, {ctci -K^c^ci +h.c.}. 
curly brackets is absent for open boundary conditions. 
Because parity of the number of fermionic excitations, 



I.e. 



of EL 



ctcfc, is a conserved quantity we can sep- 



1 '-k 



arately diagonalize hamiltonian in even and odd sub- 
sectors. In the following we will assume n is even. 
The only difference between odd and even subsectors 
is in the sign of the term coupling the last and the 
1st fermion. Diagonalization of both can be treated 
on the same footing by writing H — hY^"^^{2cjCj — 
1) + E"=i(cjcj^]^ + jCjCj+i + h.c.) and boundary con- 
ditions Cn+i — —ci for even subsector and Cn-^-i — ci 
for odd subsector. Transforming fermionic operators Cj 
with a unitary Fourier transformation to reciprocal space 



^/4 



T^^ELiC^fc exp(i^fcj), 
- 1) + 2 cos i^k)dkdl + 



fermionic operators d^, Cj — ■ 
we obtain H = h{2d\dk 

7sin(^fc)(dj,dlj. + d-kdk). In order to satisfy bound- 
ary conditions for fermionic operators Cj the allowed val- 
ues of k are k = i,|,...,n — i in even subsector and 
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2 ' 2 ' ■ ■ ■ ' " 2 
fc = 0, 1, . . . , n — 1 in odd subsector. Finally, Bogoliubov 



transformation diagonalizes two-dimensional subspace of 
dk and d-k, resulting in if = J2k^k{2flfk - 1) with 
single-particle fermionic excitation energies 



£k 



cos ( — fc) — h 
n 



7 sm ( — k) 



(15) 



In the odd subspace two fermionic eigenenergies with no 
"— fc" partner are ek=o = h — 1 and £k=n/2 = h + 1. 
Ground state of XY model is doubly degenerate [l2| for 
7^ + /i^ = 1 and this is the case occurring for Markov 
chains. Because the spectrum is (for even n) symmetric 
with respect to 0, there are also two maximal eigenstates. 
These two correspond to two invariant states of M with 
eigenvalues A = 1. For 7^ + ft,^ = 1 eigenmode ener- 
gies simplify to efe = 1 — /icos(^fc). Eigenenergies 
Ej of XY model are now obtained by filling even/odd 
number of fermionic modes. There are two eigenen- 
ergies £'1,2 = n, coming one from even and one from 
odd subsector while the second largest eigenenergy is al- 
ways from even subsector, is nondegenerate, and equal 
to — n — 4£inin = n — 4(1 — \h\ cos {n/nj). The anal- 
ysis is similar for odd n. At the end the same formu- 
las ([T5|) can be used, one only has to take all with 
a negative sign. Largest three eigenenergies are though 
given by the same i?i,2,3 = n,n,n — 4(1 — \h\ cos (tt/tt,)). 
Finally, for open boundary conditions the three largest 
eigenenergies are given for an arbitrary n by -Ei,2,3 = 
n — l,n — l,n— 1 — 2(1 — \h\ cos {n/n)), with E3 being 
doubly degenerate. 
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